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Abstract 

The production of triple Higgs ( H + H~H °), ( H + H~h °) and pair wise charged Higgs boson (H + H~) 
is studied in the context of future linear colliders within the Two-Higgs Doublet Model (2HDM) type 
II. The aim is to compare sources of charged Higgs pair through the above processes, i.e., double and 
triple Higgs production. Cross-sections are calculated at the leading order in 2HDM type II and Minimal 
Supersymmetric Standard Model (MSSM). Several orders of magnitude (~ 10 4 ) enhancement is observed in 
2HDM compared to MSSM, while no sizable enhancement is seen in muon collider versus electron-positron 
collider. The analysis is based on a heavy charged Higgs with mass above 500 GeV. It is found that double 
charged Higgs production cross-section (being the same in 2HDM and MSSM) is few femtobarns while the 
triple Higgs production can not exceed a fraction of femtobarn within the parameter space under study. 
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I. INTRODUCTION 


The main undisputed highlight of Run 1 of the Large Hadron Collider (LHC) at CERN ELS 
is the discovery of the Standard Model (SM) Higgs boson. The measured signal strengths are 
quite in agreement with SM predictions. The mass of the Higgs signal is found to be near 125 
GeV, which not only confirms the Higgs mechanism as a right approach towards giving masses to 
the electroweak particles and gauge bosons, but also puts a question of possibility of existence of 
further Higgs bosons, as it is still not clear whether the Higgs sector is indeed minimal, containing 
only a single Higgs doublet. One of the straight forward ways to address such a question is simply 
to go beyond the SM by adding a second Higgs doublet to the field content of the model. A partic¬ 
ularly well-motivated possibility along these lines is the Minimal Supersymmetric Standard Model 
(MSSM) Q] and the general (unconstrained) Two-Higgs-Doublet Model (2HDM)[4], The general 
2HDM Higgs sector contains two CP-even neutral Higgs bosons, hP, H°, a CP-odd (pseudo-scaler) 
neutral Higgs boson, A °, and pair of charged Higgs bosons, H^, whereas, hO is the SM-like Higgs 
boson and is the candidate for the signal observed at LHC. 

The purpose of this paper is to take into account all current constraints on the type-II CP- 
conserving 2HDM parameter space and determine the allowed ranges of the triple and double Higgs 
couplings to estimate their corresponding cross-sections. This work would prepare the ground for 
collider studies. Therefore two types of linear colliders, i.e., e + e _ and /r + /r _ are compared for 
all processes. Both of these type of processes have been extensively searched during the last 
years at Tevatron, Large Electron-Positron Collider (LEP) and currently at Large Hadron Collider 
(LHC). The recent results from LHC exclude a large parameter space of the light charged Higgs, 
m H ± < 160 GeV, if BR(iL ± —> tv) = 1 and heavy charged Higgs at tan/3 > 50 [ 5 ], (J. Furthermore, 
a remarkable restriction over charged Higgs masses comes from Flavour Changing Neutral Current, 
(FCNC), radiative B-meson decay, whose branching ratio BR(b —>• sj) ~ 3 x 10 -4 is measured 
with sufficient precision that becomes sensitive to new physics. The charged Higgs contribution in 
above branching ratio increases with decreasing mjj±. This channel has been studied by BaBar 
and Belle collaborations in detail {s, 9] and the up-to-date limit excludes charged Higgs lighter 
than 480 GeV at 95% C.L. H. Therefore fhe arrays preset ,hror, g hor,f fhe paper is haseri 
on m H ± > 500GeV. 

A phenomenological study of triple Higgs production, including event study and the effect of neu¬ 
tral Higgs masses, m# and m_A, on the cross section can be found in [ll] for m H ± < 500 GeV. 

In order to choose right masses for Higgs bosons of 2HDM, one should be aware that the global 
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fit to electroweak measurements requires A p to be O (10~ 3 ) [ 12 |. This requirement does not allow 
large mass splitting between Higgs bosons. Therefore, we adopt degenerate masses for neutral 
Higgs bosons, i.e., m# = tha = m H ±. 

Regarding the neutral Higgs mass recently CMS and ATLAS experiments have excluded a 
wide range of m H o in MSSM via a study of H —>• rr channel at yfs = 8 TeV 3 , Rjj]. Higgs boson 
decays to gauge bosons such as H —>• WW, ZZ and 77 have not been studied for MSSM neutral 
Higgs bosons. They are, however, considered for the light SM Higgs boson as different sources of 
Higgs boson production. Our chosen MSSM points are outside the LHC excluded area if one sets 
mu = mA = fn H ± > 500 GeV and tan/3 = 10. 

This study is intended to be suitable at e + e _ colliders or muon colliders. The muon collider is 
expected to get the integrated luminosity around 125 fb~ 1 at yfs = 1.5 TeV and 440 fb~ 1 at yfs 
= 3 TeV. This is a unique machine, to be designed to provide high luminosity, very small energy 
spread, excellent stability, and good shielding of muon beam decay backgrounds. 


II. THEORETICAL FRAMEWORK 


The theoretical basis is a two Higgs doublet model with the general potential as follows 


0,Q. 


V = + m\ 2 ^\^2 — 771^2 + h.C. 


+ ^Ai + 1 a 2 

+ A 3 (V^l) ) + A 4 (^1^2) ( 




+ -A 5 ( 4 >{$ 2 ) + A 6 + A 7 


$i$ 2 ) + h.c. 


The free parameters of such a model are 


Ai, A 2 , A3, A4, A5, A6, A7, m\ 2 1 tan P 


(1) 


in the general basis. The CP violation or Flavor Changing Neutral Currents (FCNC) are not 
assumed as they are naturally suppressed via the Natural Flavor Conservation (NFC) mechanism 
by imposing a Z 2 symmetry on the Lagrangian [17, 18] leads to the A 6 = A 7 = 0. Furthermore we 
choose sin(/3 — a) = 1 which takes the region of study very close to MSSM parameter space. With 
the above setting, i.e., tan (5 = 10 and sin(/3 — a) = 1, free parameters can alternatively be taken 
as: mh,mH,mA,mH±,mi 2 ,tan/3 because there is a correspondence between Higgs boson masses 
and A values. We use 2HDMC package [19] to ensure that chosen parameters are consistent with 
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FIG. 1: Higgs boson masses and the value of m \2 which respect physical requirements on model potential 
as well as experimental limits and observations on Higgs boson masses. 

current experimental limits and respect also the potential unitarity, perturbativity and stability. 
A point is chosen if it satisfies all above requirements. Figures [T] and [2] show Higgs boson masses 
used in the analysis and their corresponding A and m .12 values extracted from 2HDMC. The choice 
of the default set of MSSM parameters are Msusy = 1000 GeV, // = 200 GeV, A t = 1000 GeV, 
A/, = 1000 GeV and A r = 1000 GeV, which reflect the benchmark m™ ax scenario. 


III. CHARGED HIGGS PAIR PRODUCTION 


Double Higgs or pair production of charged Higgs bosons at linear colliders, has been extensive 
investigated in the literature mainly in MSSM in the context of e + e _ and /r + [i~ colliders, 20 2c 


The two types of charged Higgs decay have been adopted in those analyses, i.e., H ± — >• tv and 
—>• tb. One of the week points of this channel is that it is limited by the center of mass energy of 
the collider and a charged Higgs mass heavier than y/s/2 is not produced unless a small negligible 
rate due to off-shell production is considered. In Figure[3]the total production cross-section cr(/x/r —» 


4 







FIG. 2: The values of A parameters which result in Higgs boson masses shown in Fig. [T] 

H + H~ ) (in pb) is plotted as a function of center-of-mass energy yfs (in GeV) for degenerate masses 
of Higgs bosons in 2HDM type II. The charged Higgs pair production cross section is independent 
of tan (3 and tua values. As is seen in the Figure [3j the corresponding production rates reach few 
femtobarns, attaining several thousand events per 1000 fb 1 of integrated luminosity. It should be 
noted that there are Feynman diagrams which could result in enhancement of the production rate. 
They are s— channel diagrams containing neutral Higgs bosons contributions and their couplings 
with leptons as well as t —channel charged Higgs diagrams shown in Fig. SI However such diagrams 
have very small contributions to the total cross section and the Z/ 7 —mediating diagrams dominate. 
Since the relevant coupling in this diagram is ZH + H~ (e. cot 2/3) and r )H + H^ (e) and both of them 
are model independent gauge couplings, the obtained cross sections are the same in 2HDM and 
MSSM. Table Q] compares contribution of different diagrams in the total production cross section 
at a certain point in parameter space. The conclusion for this observation is that double charged 
Higgs production (H + H~) is produced at /x + /z - and e + e _ colliders with no sizable difference. 
Since production cross section is the same at both 2HDM and MSSM, it is difficult to distinguish 
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FIG. 3: Total cross-section a(H + H ) ( in fb) as a function of \fs for the tree-level Higgs boson pair 
production at 2HDM or MSSM at a muon collider. 



FIG. 4: The possible Feynman diagrams for double Higgs production at lepton colliders. 



H+ 

MSSM 

H~ 

2HDM 

Collider type 

Zh 

h°/H° 

Total 

Z/l 

h°/H° 

Total 


2.73 

1.40 x 10" 6 

2.73 

2.74 

1.41 x 10 -6 

2.74 

e + e~ 

2.73 

1.27 x 10~ 13 

2.73 

2.74 

7.61 x 10" 16 

2.74 


TABLE I: Total cross section (in fb) of double charged Higgs at e + e and colliders at yfs = 3 TeV. 

The Higgs boson masses are run = tua = rriH± = 500 GeV. 
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the two models using double Higgs ( H + H ) channel. It should be noted that one can look to 
single charged Higgs production in association with W boson. This channel has been studied in 


24-l2g] and shows that signal cross section is larger at muon colliders and is almost independent 


of the charged Higgs mass and increases with tan/3 . A similar analysis of this channel at LHC 
has been repoted in [27]. In the rest of the paper, triple Higgs production is studied. It might 
be possible that the triple Higgs boson self-interactions play a role in this endeavor, however, as 
is shown in the next sections, cross section of triple Higgs production [H + H~H °, H + H~h°) is 
smaller than double Higgs production ( H + H ~). 

IV. TRIPLE HIGGS PRODUCTION 


The trip 


e Higgs production has been studied in different papers in the context of linear 


28-[32]. The production processes under consideration are H + H H° and H + H h° 


colliders 

in fj, + fx~ annihilations within the 2HDM type II with their corresponding trilinear Higgs 
bosons couplings given in the Eqns. [2] and [3] together with their corresponding MSSM values 
in Eqns. [4] and [5l where usual abreviations such as sw = sin^vv and sp = sin/3 have been used [33]. 


H ± H ± H°(2HDM) : 


2 ( cis 2 ps a m 2 H + c a s 2/ 3 simj 1 - 2s Q+ p[i 2 12 + cp- a s 2 20 m 2 H ±) (2) 

mw.sw-s 2/3 


H ± H ± h°(2HDM) : 


---—(2 c a+0 fi( 2 - c a CpS 2 pml + s 2 ps a spm 2 h - s 20 sp- a m 2 H ±) (3) 

.Syj/ *^2/3 


^^H^iMSSM) = C ' mz 

2mws\v 


(C2W C/3—a + S 2 ps a+ p ) 


( 4 ) 


H^H^h^MSSM) = (2 c 2 w sp- a +c 2 ps a+ p) 


—e.m 


( 5 ) 


2m,]ysw 

Figure [5] shows Feynman diagrams of triple Higgs production. In this case there is tan /3 dependence, 
contrary to the case of H + H~ production. The H + H~H° coupling effectively grows as tan/3 or 
cot/3 for tan/3 >> 1 or tan/3 <<1 respectively. The corresponding cross-section can vary either 
by tan 2 /3 at larger tan/3 values or by cot 2 /3 at small tan/3 values respectively in contrast to the 
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situation in MSSM, where the triple Higgs coupling undergoes radiative corrections m and does 
not have any possible source of enhancement as can be seen from Eqn|3] showing tree level coupling 
in MSSM. The Eqn. [I] clearly indicate that the couplings are naturally gauge like and hence the 

2l|]. Figures [6] and [7] show cross sections of H + H~H° 


expected cross-section remains rather small 


and H + H~h° production for different set of charged Higgs masses and center of mass energies. 
One can see the contribution of different parts of the set of diagrams and compare them in a table 
as was done for double Higgs production. Table [TT] shows such information for H + H~H° and 
H+H~h° processes, where a comprehensive comparison between 2HDM and MSSM results, for 
/r + /x - and e + e _ colliders has been done. 

At this stage, a comparison between the two types of processes, i.e., double and triple Higgs 
production, is performed by finding a center of mass energy which gives the maximum allowed 
cross-section. In Table HTT1 the cr max (H + H~H°), (7 max ( H + H~h° ) and o(H + H~) at two different 
center of mass energies (yfs = 1 TeV and yfs = 1.5 TeV) are listed. The abreviated name ”NP” 
stands for ’’Not Possible”, due to not enough center of mass energy to produce a given event. As 
is seen from Table m the MSSM cross-sections for both triple Higgs cases are extremely small 
reaching the largest value of a max (H + H~ H°) « 10 -6 fb and a max {H + H~h°) ~ 10 -5 fb. 
Comparing 2HDM and MSSM cross sections, it is observed that one could reach a factor of 10 5 
enhancement in 2HDM cross section of H + H° , while for H + H~h°, the 2HDM cross section is 
typically 10 times larger than the corresponding MSSM value. On the basis of these results, the 
triple Higgs boson channels are generally much more promising in 2HDM than MSSM framework, 
but they do not reach the level of H + H~ cross section. The double Higgs production cross-section 
is substantially larger than triple Higgs channels either in MSSM or 2HDM. Therefore any study 
of the charged Higgs at linear colliders should be based on double Higgs production, unless a large 
integrated luminosity and center of mass energy is available for production of triple Higgs processes. 
Such events can only be considered as complementary processes to shade some light on triple Higgs 
couplings. 


V. CROSS SECTION AT e+e' AND //+/C COLLIDERS 

In this section, cross sections of three selected channels at two possible choices of colliders, i.e., 
e + e _ and are compared. As shown in Figure 01 the charged Higgs pair production may 

proceed through three types of diagrams. The left diagram includes s-channel electroweak propa¬ 
gators while the middle one consists of neutral Higgs bosons as propagator. The CP odd neutral 
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FIG. 5: The possible Feynman diagrams for triple Higgs production at lepton colliders. 
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FIG. 6: Triple Higgs production ( H + H H°) within 2HDM type II for various m#± values. Neutral Higgs 
bosons have masses equal to the mass of the charged Higgs. 


9 














0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0 


FIG. 7: Triple Higgs production ( H + H h°) within 2HDM type II for various m#± values. Neutral Higgs 
bosons have masses equal to the mass of the charged Higgs. 



FIG. 8: Branching ratio of charged Higgs decays in 2HDM type II. 
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H+H~H° 

H+H~h° 

Collider type 

Model type 

zh 

h°/H° 

Total 

zh 

h°/H° 

Total 


MSSM 

1.93 x 1CT 7 

3.49 x 10~ 6 

3.69 x 10 -6 

1.69 x 10" 5 

2.85 x 10~ 8 

1.69 x 10" 5 


2HDM 

0.116 

2.24 x 10“ 6 

0.116 

4.44 x 10 -4 

1.90 x 10 -10 

4.44 x 10“ 4 

e + e~ 

MSSM 

1.93 x 10" 7 

7.77 x 10- 11 

1.93 x 10" 7 

1.69 x 10" 5 

6.34 x 10~ 13 

1.69 x KT 5 

e + e~ 

2HDM 

0.116 

5.10 x 10" 11 

0.116 

4.44 x 10“ 4 

1.03 x 10" 19 

4.44 x 10" 4 


TABLE II: Total cross-section (in fb) comparison between electron-positron and muon-muon colliders as 
well as among MSSM and 2HDM at \fs = 1.5 TeY. 


Process 

m. H ± 

CTmax (1.5 TeV) 

Umax (3 TeV) 

<?max (1.5 TeV) 

<?max (3 TeV) 


GeV 

MSSM 

MSSM 

2HDM 

2HDM 

H + H~ H+H~ 

500 

5.26 

2.72 

5.26 

2.72 

fji+fj,- H+H~H° 

500 

NP 

3.68 x 10" 6 

NP 

0.12 

H+H~h° 

500 

1.44 x 10" 5 

1.69 x 10 -5 

3.73 x 10- 4 

4.44 x 10" 4 

H+yT H+H~ 

700 

5.28 

2.25 

5.28 

2.25 

Ai+M“ -> H+H~H° 

700 

NP 

3.86 x 10 -7 

NP 

1.81 x 10" 2 

H+V~ H+H~h° 

700 

NP 

1.13 x 10 -5 

NP 

2.2 x 10~ 4 

IL + H~ H+H~ 

900 

NP 

1.66 

NP 

1.66 

fj+H~ H+H~H° 

900 

NP 

3.86 x 10~ 8 

NP 

5.77 x 10" 4 

H+liT H+H~h° 

900 

NP 

5.27 x 10 -6 

NP 

1.02 x 10" 4 


TABLE III: The maximum cross-section (in fb) for the leading order double Higgs and triple Higgs processes 
within MSSM and 2HDM at two different center of mass energies y/s =1.5 TeV and yfs = 3 TeV. 

Higgs boson does not contribute in the production cross section due to having zero coupling with 
charged Higgs pair. Therefore we only take into account the electroweak and CP-even neutral 
Higgs bosons contribution. The amount of cross section coming from each diagram may be con¬ 
siderably important probably after a long run of any one of these two colliders in future. There is 
also a t-channel process (the right diagram in Fig. [4]), which has negligibly small contribution at 
both e + e and n + /j, colliders due to the small Yukawa coupling between the charged Higgs and the 
leptons. The Table U demonstrates that H + H~ cross section is dominated by the Z /7 mediating 
diagram (the left diagram in Fig. [4]) and no difference is obtained between e + e” and even 

though by introducing the non zero masses of electron and muon and their corresponding couplings 
in the Lagrangian. There is in fact a large enhancement in h/H mediating diagram, however, the 
total cross section is not affected by this enhancement due to the negligible contribution of collider 
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dependent diagrams. 

The same conclusion holds for the triple Higgs production processes shown in Figure EJ The di¬ 
agram with only Z /7 propagator (Fig. [5]-b) has the dominant contribution and as seen in Table 
m the effect of collider dependent diagrams, i.e., the ones with neutral or charged Higgs coupling 
with leptons, is small. 

In summarizing the results we can conclude that for the charged Higgs searches the most suitable 
channel from cross section point of view is charged Higgs pair production having largest cross- 
section of the order of few fb in MSSM and 2HDM, while the selected triple Higgs processes 
(H + H~H °, H + H~h° ) have very tiny cross-sections. The triple Higgs searches could benefit from 
the neutral Higgs decay to bb which can be identified by 6 —taging algorithm. However, due to 
the larger particle multiplicity compared to H + H~ production, selection efficiencies of the triple 
Higgs are expected to be small. The final state to analyze for both double and triple Higgs cases 
could be defined by selecting one of the two possible decay channels of the charged Higgs, i.e., 
H± tv and H± —» tb. As Figure [ 8 ] shows, the charged Higgs branching ratio of decay to tv goes 
down to the level of 0.1, for m#± > 500 GeV, while decay to tb receives a branching ratio of 0.9. 
Since cross sections are small, especially, for triple Higgs production, the main channel to study is 
expected to be the one with top and bottom quarks in the hard scattering final state. Of course 
in this case a large number of particles is produced due to the decay of the top quark to W boson 
and its subsequent decay to leptons or light jets. One the contrary, the final state with r leptons 
involved, may be identified using a r identification algorithm but it certainly suffers from the very 
tiny number of events in hand. In any case the double Higgs production is expected to be more 
promising than the triple Higgs production due to the larger cross section. 

VI. SUMMARY 

The most promising channels were studied for the charged Higgs production through double 
or triple Higgs production at future e + e _ and linear colliders in the framework of CP- 

conserving Two Higgs Doublet Model type II and Minimal Supersymmetric Standard Model. The 
sizes of double and triple Higgs production cross-sections were computed both in MSSM and 2HDM 
by respecting all the current experimental and phenomenological constraints on Higgs masses, re¬ 
quirement of perturbativity and vacuum stability condition. The maximum cross-sections at two 
different center of mass energies of 1.5 TeV and 3 TeV were obtained to compare these processes 
between MSSM and 2HDM where several order of magnitude difference was observed. A com- 
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parison between e + e~ and fj. + /x - colliders was also performed leading to no sizable difference at 
both colliders. The conclusion is that double Higgs production cross section is much larger than 
triple Higgs, thus keeping it as the main sourse of the charged Higgs at linear colliders. There 
are diagrams which depend on the collider type (incoming lepton) but they do not have a sizable 
contribution to affect the total cross section. This fact is observed for both double and triple Higgs 
productions. The double Higgs production is less limited by the center of mass energy of the system 
compared to the triple Higgs production and thus can probe heavier charged Higgs bosons at a 
given center of mass energy. The final state to choose for a data analysis would be H + H~ —>• tbtb 
and H + H~H°(H + H~h°) —> tbtbbb for the top and bottom final states, and, H + H~ —>• T + VTU 
and H + H~H°[H + H~h°) —>• T + vTvbb for t and missing transverse energy final states. While final 
states with t and missing transverse energy are suitable for light charged Higgs searches, a heavy 
charged Higgs with mjj± > 500 GeV would require a final state including top and bottom quarks. 
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